Warped Universal Extra Dimensions 



Am'bal D. Medina and Eduardo Ponton ^ 

"'Department of Physics, University of California, One Shields Ave. Davis, CA 95616, USA 
^Department of Physics, Columbia University, 538 W. 120th St, New York, NY 10027, USA 

June 3, 2011 

Abstract 

We consider a 5D warped scenario with a KK-parity symmetry, where the non- 
trivial warping arises from the dynamics that stabihzes the size of the extra dimension. 
Generically, the lightest Kaluza-Klcin (KK) particle is the first excitation of the radion 
field, while the next-to-lightest Kaluza-Klcin particle is either the first excitation of 
the (RH) top quark or the first KK-parity odd Higgs. All these masses are expected 
to be of order the electroweak scale. We present simple analytical expressions for 
the masses and wavcfunctions of the lowest lying KK modes, and derive the Fcynman 
rules necessary for phenomenological applications. The framework allows to interpolate 
between a strongly warped scenario a la Randall- Sundrum (RS), and a weakly warped 
scenario that shares properties of both RS and Universal Extra Dimensions models. 



1 Introduction 



Warped extra dimensional scenarios provide an appealing solution to the hierarchy prob- 
lem In addition, when the SM gauge bosons pi [3] and fermions [1] arise from higher- 
dimensional fields one gains a rather interesting understanding of the SM fermion hierarchies 
and CKM mixing angles, based on the localization of these modes along the extra dimen- 
sion [5] : light fermion fields are localized near the "UV brane" while heavier fermions, most 
notably the top quark, are localized closer to the "IR brane". At the same time danger- 
ous fiavor-changing neutral current effects are suppressed, and can be roughly consistent 
with low-energy constraints even when the new physics is at the TeV scale. The associated 
spectrum of Kaluza-Klein (KK) modes provides exciting signals at the LHC. 

It was recognized early on |6l |7] that a complete solution to the hierarchy problem requires 
the specification of dynamics that stabilizes the size of the extra dimensions. In the Randall- 
Sundrum (RS) scenario [I], this requires that the interbrane separation be stabilized at a 
factor of ~ 35 in units of the inverse curvature scale. This can be achieved either at tree-level, 
by considering for example a SM singlet scalar that acquires a non-trivial vacuum expectation 
value (VEV) [8], or at the quantum level (see e.g. [TOl HH [El [131 [H]). An important 
consequence is the presence of a 4D scalar mode, called the "radion" , that corresponds to the 
fiuctuations in the interbrane distance. The radion mode is expected to be parametrically 
lighter than the KK states, and therefore potentially important in the phenomenology of 
these scenarios. The generic couplings of the radion to SM fields were worked out in Refs. [151 
[T6l [T71 [18] . Except for a possible mixing term with the Higgs field [HI [201 HI], these 
interactions are non-renormalizable and suppressed by a "radion decay constant", A^. 

In this paper we propose stabilizing the size of the extra-dimension, parameterized by 
y G [—L,L], in such a way that the non-trivial warping arises from the backreaction of a 
Goldberger-Wise (GW) scalar, with the geometry being symmetric about y = 0. As a result, 
the boundaries of the space, at y = ±L, are IR boundaries were the warp factor is smallest, 
while the maximum of the warp factor is achieved at y = 0. Thus, although no UV brane is 
put in by hand, a dynamical "UV brane" is generated in the middle of the extra dimension. 
In fact, in a certain limit this UV brane becomes infinitely thin, and the framework reduces 
to two copies of the RS scenario P glued together at the UV brane, as considered in [22j. 
However, generically, the dynamical UV brane is "fat", and the geometry is not simply AdSs. 

One important consequence of the setup, that goes beyond the most studied implementa- 
tions of the RS proposal (for a review, see e.g. pS]), is the existence of a discrete Z2 symmetry 
under which all the even-level KK states (which include the zero-modes) are 2'2-even, while 
the odd-level modes are Z2-odd (i.e. a KK parity similar to the one present in Universal 
Extra Dimensional scenarios or UED's [2^ [25]). As a result, the lightest KK-parity odd 
particle (the LKP) is stable, and a potential dark matter (DM) candidate (for an alterna- 
tive proposal, see [2S])- However, unlike in UED scenarios [23 [2H], the KK states are not 
necessarily highly degenerate, with the degeneracy broken by loop-level effects, but instead 
larger mass splittings can be expected (although there are some interesting degeneracies, 
related to localization effects). A similar idea of a warped framework with a KK parity was 
explored in [22] , where the viability of the first KK excitation of the Z gauge boson as a DM 
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candidate was investigated. Here, we point out that the LKP is generically expected to be 
the first KK excitation of the radion field (more precisely of the radion/GW-scalar system). 
We show that its mass is highly degenerate with the radion mass, and that its interactions 
are also controlled by A^. Then, the question of the viability of the first KK radion as a DM 
candidate naturally arises, and answering it becomes essential in such warped frameworks 
with a KK parity. The point is that, although it is expected that the LKP has a mass of 
order the EW scale, its interactions may not be sufficiently strong to ensure the annihila- 
tion of KK radions during the early history of the universe to levels below the "overclosure 
limit". The WIMP miracle is not necessarily automatic, and the danger of overproducing 
the LKP's becomes a pressing issue. Fortunately, there are a number of scenarios one can 
envision where the KK radion relic density can be suppressed and explain the observed DM 
density. We explore these in detail in the companion paper [29] . 

In this work we lay out the formalism necessary to analyze the properties of the LKP, 
which includes also a discussion of the KK excitations of fermion and gauge fields. Since we 
consider a rather general gravitational background, it is difficult to find closed solutions for 
the KK wavef unctions. However, we obtain reasonably accurate analytical approximation 
for the lightest states in the theory, which include the radion, the KK radion and possibly 
the first KK excitation of the (RH) top. These become more accurate in the limit that 
the UV brane is fat, which corresponds to an interesting deformation of the RS proposal. 
We also point out that fermion zero-mode localization is naturally achieved by coupling the 
5D fermions to the GW scalar that induces the warping and stabilizes the size of the extra 
dimension. Thus, no ad-hoc 5D fermion masses need to be introduced. Regarding the EWSB 
sector, it is shown that it is naturally a two-Higgs doublet model (THDM), with the second 
doublet being a KK-parity odd excitation of the SM Higgs. 

Our formalism also allows us to interpolate between a strongly warped limit (a la RS) and 
a limit where the curvature is small, which is similar to the UED assumption. An interesting 
intermediate case arises when the curvature is small compared to the 5D Planck mass, but not 
exactly zero. In this case, the hierarchy between the Planck and the electroweak (EW) scales 
arises in part from a moderate warp factor, and in part from the smallness of the curvature 
scale. Such a scenario can share phenomenological properties of both RS and UED scenarios. 
It can also give rise to KK radion DM as a non-thermal relic |29] . Our approach is somewhat 
different from the recently studied "soft-wall" scenarios [30]-[39], but it could be interesting 
to combine the two ideas to obtain a fully dynamical warped compactification. 

The outline of this work is as follows: in Section |2] we give several examples that illustrate 
the stabilization of the extra dimension with an automatic KK-parity symmetry. We explore 
various limits, imposing a minimal number of phenomenological constraints, and define two 
benchmark scenarios to be used as reference. In Section [3} we discuss the KK decomposition 
of the radion/GW-scalar system, thus defining the radion decay constant A^. We also give 
simple analytical expressions for the mass and wavefunction profiles of both the radion 
and its first KK excitation (the LKP). Section |4] is devoted to a study of fermions in the 
Z2 symmetric backgrounds. We point out that in the fat UV brane limit it is possible 
to characterize the localization of the zero-mode via a dimensionless c-parameter that is 
analogous to the treatment in the pure AdSs background, although there are important 
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differences. We also give approximate analytical expressions for the mass and wavefunction 
of the first fermion KK resonance in cases where it is parametrically lighter than the KK scale 
(hence the likely next-to-lightest KK-parity odd particle, or NLKP). Section [s] is devoted to 
bulk scalars as applied to the Higgs field. We also present the Higgs Yukawa interactions 
with bulk fermions. Gauge fields are briefly discussed in Section [6] In preparation to the 
DM analysis to be presented in [29] , we work out in Section [T] the Feynman rules for the 
interactions between the radion/KK-radion and the KK fermions, gauge bosons and the 
Higgs doublets. We summarize and conclude in Section |8| 



2 General Setup 

In this section we describe the stabilization mechanism that leads to a non-trivial Z^- 
symmetric warping. We illustrate this by choosing scalar potentials that allow for closed 
analytical solutions that take into account the scalar backreaction on the metric, but em- 
phasize that the properties we are interested in do not depend on any given particular choice 
of potential. We discuss an initial set of phenomenological constraints, and define two bench- 
mark scenarios, one with strong warping and one where the curvature is small compared to 
the 5D Planck scale. 



2.1 Backgrounds with Z2 Symmetry 

We start with a 5D real scalar $ minimally coupled to gravity according to 



S 



M 



dM 



where M5 is the (reduced) 5D Planck mass, TZ^ is the 5D Ricci scalar, and the last term 
allows for operators localized on the boundary of the fifth dimension (which will be specified 
in the next subsection). We assume that the fifth dimension is compactified to an interval, 
parameterized by ?/ G [—L,L]. Our interest is in dynamics that stabilizes the size of the 
extra dimensions, 2L, while leading to a background that is symmetric about y = 0, so that 
a geometric Z2 symmetrjj^is present, as in [22]. However, unlike Ref. [22] we do not have a 
central brane at y = 0, although as we will see, in certain limits, such a brane arises from 
the scalar dynamics. 

We are interested in solving for the coupled gravity /scalar system, taking into account 
the backreaction of a scalar VEV on the geometry, since this results in the stabilization of 
the extra dimension a la Goldberger-Wise [H]. We restrict to backgrounds exhibiting 4D 
Lorentz symmetry, which can always be put in the form 



ds' 



'^^^y^^.udx^dx" - dy'' , 



(2) 



^As established in [30], Chern-Simons terms do not break the KK-parity symmetry, so that the questions 
raised in [41 for httle Higgs scenarios with T-parity do not apply. 
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where y is the proper distance. The weak energy condition imphes that A"{y) > |12] so 
that e~^^y^ must be a convex function and the two boundaries at y = ±L are IR branes 
in the RS terminology. As shown in Refs. ^\ 1^ . Einstein's equations, together with the 
scalar equation of motion, can be solved exactly for the class of scalar potentials given by 



where is an arbitrary "superpotential" [W has mass dimension 4, while $ has mass 

dimension 3/2 in 5D). For the present purpose, the important point is that the scalar back- 
ground, (f){y) = {^{y)), can be obtained by solving the first order differential equation]^ 

and that once the solution (f){y) is found, A'{y) is simply given by 

A'{y) = ^i^i^iy))^ (5) 

which can be integrated immediately to yield y4(?/)|^ Eq. Q shows that (f)' is even under 
y — )• —y, provided dW/d(f) is even -hence W{(j)) is odd- under — )■ —0. It then follows from 
Eq. ([5]) that A'{y) is odd, and therefore A(y) is even under y — )■ —y, as we want. Thus, 
imposing the Z2 reflection symmetry allows us to focus on odd superpotentials. We also 
see that for the stabilizing scalar, $, the discrete symmetry should be a combination of the 
geometric reflection y —y and an internal parity symmetry: the action of the relevant 
Z2 on the GW scalar is — t- —^(x^, ~y)- The scalar potential derived from Eq. ^ 
is clearly even under this Z2, and the latter is respected by odd "kink-like" profiles for the 
VEV, 0(1/). 

For instance, the simplest ansatz for the superpotential is a linear one: 

Wi{^) = 2m0o$ , (6) 

where, for later convenience, we have parameterized the coefficient in terms of a mass scale 
m and a scalar "VEV" 0o (with mass dimension 3/2). Both Eqs. ^ and ([s]) can be trivially 
integrated, and give 

(f){y) = (j)omy , (7) 



A(y) - 



2 



■^Eqs. Q and ^ guarantee that the equations of motion for the gravity/scalar system are satisfied when 
the potential is given by Eq. ([s]); see [U for details. 

•^Alternatively, given A{y) one can combine Eqs. (jij) and ^ using = {l/4>')dy to obtain (j)' = 
±-y/f 2M|A"(?/), which can be integrated to yield 4>{y), up to a sign. If 4>{y) can be inverted to obtain 
y = y{(j)), the superpotential is given by W{(j)) — 6M| A' (?/((/))). 
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where the condition (j){—y) = —(piy) sets the integration constant to zero in Eq. ([T]), and 
we choose ^(0) = by rescahng the 4D coordinates by the appropriate constant scale 
factor. This simple example will turn out to capture the correct physics of a large class of 
superpotentials, once the radion is stabilized and a hierarchy is induced. 
A richer example arises from a cubic superpotential: 

= ^(<^o$-^^') ' (9) 

with m an arbitrary mass scale that, without loss of generality, can be chosen to be positive. 
In this case, Eq. Q reads 



= - [0^ - <t>iyr] , (10) 

which can be recognized to give the standard fiat space "kink" solution |^ 

(j){y) = (j)otanh{my) , (11) 

with a thickness of order 1/m and asymptotic VEV ±0o- Following the above described 
procedure, one then finds 

k 

^(y) = ^ — [tanh^(m?/) + 41ogcosh(my)] , (12) 

where we defined the scale A; by 0g = 9A;M|/(2m), and again we chose A{0) = 0. Unlike the 
example defined by the linear superpotential ([6]), this model has the virtue that for \y\ ^ 1/m 
one has A{y) — > k\y\ + const., so that asymptotically one gets AdSs with curvature given by 



k. In fact, it is clear that this corresponds to the limit where the domain-wall (11) is very 
narrow, i.e. m ^ k, and models a UV brane at y = 0, so that the RS solution is recovered. 
We note that the potential that follows from the cubic superpotential (|9|) is 



9-n of 3k\ o 4:'m'^ f m\ . 4m^ 



y ($) = -mkM^ - 1 + — + —p. 1 + ^ $4 ■ (13) 

^ ^ 4 ^ \ mj 9M| V AkJ 2A3kMl ^ ' 

From an effective field theory perspective, this is a scalar potential subject to a parity 
symmetry $—)■—$, and truncated at sixth order in $. However, it is non-generic in two 
ways. First, the constant term has a fixed value, which is necessary for the ansatz of Eq. ^ 
-with flat 4D sections- to hold. In addition to this well-known Cosmological Constant 
Problem, we see that the potential contains three operators that are controlled by only two 
parameters, thus seemingly implying an additional fine-tuned relation (at least in the absence 
of supersymmetry) . The two parameters can be chosen as k/M^ and m/M^ (M5 simply sets 
the overall scale, and is associated to the strength of the gravitational interactions). However, 



4r 



The scalar potential, Eq. (13 1, depends only on 0q. The sign of parameterizes the two solutions 



related by y — > — y, which correspond to a positive/negative VEV ai y — +od. 
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we stress that this additional relation plays no essential role in our applications. Indeed, we 
see that in the thin domain-wall limit, m 3> fc, and using < (po, the ratio of the 

$^ operator to the operator in Eq. (13) is of order k/m <^ 1. Hence, the dynamics is 



controlled by the quadratic and quartic operators, the last term in Eq. ( 13 ) representing only 
a small perturbation. In the opposite limit, m — )■ 0, and for |0| ~ 0o, the three operators are 
parametrically equally important. Nevertheless, in the central region, where \(f){y)\ ^ (po, 
the dynamics is actually controlled only by the quadratic term, thus reducing precisely to 
the case of the linear superpotential of Eq. ([6]). As we will see, this will be important in the 
following and makes our results rather generic. The point is that we are using the special 
class of potentials that can be obtained from a superpotential via Eq. ^ only to obtain 
simple, closed solutions that fully take into account the backreaction of the scalar on the 
geometry, but our physical results will not depend on this simplifying assumption. 

Our final example is based on the observation that the cubic superpotential above corre- 
sponds to the first two terms in the Taylor expansion of a sine: 



V2m(f) 

n sin 



V2 



(14) 



This leads to a scalar profile 



0(l/) = y/2(j)Q tan ^ 



tanh 



my 



.V2j\ ' 

that is also a domain wall of width of order 1/m, centered at ?/ = 0. Correspondingly, 

k 



(15) 



m 



log 



cosh ( \/2my 



(16) 



where now 0o and k are related by 0o = 6A;M|/(\/2m), and we again made sure that the 
warp factor is normalized to unity a.t y = 0. As in the second example, the geometry 
is asymptotically AdSs with curvature k, i.e. one has A{y) — )■ k\y\ + const, for \y\ ^ 
1/m. Again, the IR-UV-IR model of Ref. [22] arises in the limit m ^ k. Our comments 
above also apply to this example: although the potential that follows from Eq. ^ with the 
superpotential (14) contains an infinite tower of operators, only a few of them are actually 



relevant in determining the background solution: the quadratic and quartic terms in the "UV- 
brane" limit, m 3> A;, and the quadratic term in the central region y < 1/m. Nevertheless, 
the "sine superpotential" has the advantage that it allows us to interpolate between narrow 
and wide domain-walls with a simpler warp factor than for the cubic superpotential. For 



concreteness, in the numerical studies of the rest of this work we will use Eqs. (15) and (16) 



with the understanding that the results differ from the other (and similar) examples only in 
minute details. 



^The asymptotic VEV's are ±(7r/2\/2) w ±1.110o- 
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2.2 Radion Stabilization 



We described in the previous subsection the bulk solutions that follow from a given superpo- 
tential. However, one must ensure that these solutions are actually allowed by the boundary 
conditions at y = ±L. One can see immediately from Eqs. Q and ^ that this requires 



1 dW{(j)) 

2 86 



A'i±L) 



y=±L 



6M| 



(17) 



y=±L 



On the other hand, the set of consistent boundary conditions can be found by varying the 
action Eq. ([T]) with respect to both $ and A, and requiring that the surface terms obtained 
from the integrations by parts vanish. Assuming that £4 in Eq. ([T]) does not involve y 
derivatives, e.g. if it is a pure potential term, this leads to [H] 



<I>'(±L) 



1 dCi 

2 ~d¥ 



A'i±L) 



y=±L 



6M| 



y=±L 



Following Refs. [HI |T6] one can enforce that the scalar field attain a given value 0o (iii 
absolute magnitude) at the boundaries by writing an even boundary term |^ 

A£4($) = -7($'-0^)' , (19) 



that does not contribute to Eqs. (18) in the limit 7 — > +00. Thus, we see that one can use 

W^'(0o)($-0o) + A£4($) , 



A($) 



-W^(0o 



(20) 



which leads to Eqs. (17) [for 7 — +00], as required. In the limit of a fixed VEV, 0o, at the 
boundaries, we must have 



0(L) = y20otan"^ 



tanh 



mL 

71 



(21) 



where the second equality holds for the case of the "sine" superpotential, Eq. (14). This 
equation fixes the size of the extra dimension to be 



mL = \/2tanh 



tan 



a/200 



and depends only on the ratio (f)o/4>o- Then Eq. (16) gives 

k 



A{L) 



m 



log 



cosh ( \/2mL 



(22) 



(23) 



^One can also write AC4 = —7 ($ — (f>o) at y ~ +L, and = —7 ($ + 0o) at j/ = — L, so as to 

preserve the $(j/) — > —^(—y) symmetry. In this case, the constant 7 has positive mass dimension and 
perhaps it is more natural than in the case of Eq. ( 19 ) for it to take a large value. It also has the virtue 



that it selects the "domain-wall scalar profile" over a constant profile. In any case, the most important 
(simplifying) assumption is that the value of the scalar field be frozen on the boundaries. 

^Notice that if ^o/^o is such that | tan ^ J^^ ^ | > 1 there are no (real) solutions for L from Eq. (22 1. 

This is simply the statement that the brane VEV, (f>o, cannot be larger than the asymptotic VEV of the 



bulk solution, Eq. (15), for the configuration to be allowed. 
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which is parametrically large for m <^ k, i.e. away from the "thin UV brane" hmit, even 
when mL ~ 0{1). It is also large in the limit that the extra dimension is large, mL 3> 1 
and kL ^ 1, although this requires the brane and bulk VEV's to satisfy 0o/("\/20o) ~ 7r/4, 
a tuned limit. 



2.3 From Large to Small Warping 

The model described above depends on three dimensionless ratios: k/m, 0o/0o and k/M^. 
The latter ratio determines the degree of warping (from nearly flat when k <^ M5 , to strongly 
warped when k ~ M5). The case of small warping may require fine-tuning to make the 5D 
cosmological constant much smaller than the fundamental scale M5 (in addition to the fine- 
tuning associated with the vanishing of the 4D CC). Nevertheless, we will consider such a 
small warping limit since it allows us to make the connection with the well-studied models of 
Universal Extra Dimensions (UED) [24j, which assume that the 5D curvature is negligible, 
while still showing some interesting differences when k is small but non-vanishing. In the 
opposite limit of strong warping we can place a rough upper bound on the size of the curvature 
from the requirement that there be some modest hierarchy between k and the cutoff of the 
theory, A5, so that neglecting higher- dimension operators in the Einstein-Hilbert action is 
justified. Based on NDA in extra-dimensional theories [IS] we expect that ~ MI//5, with 
^5 = 24:71^, a 5D loop factor. This suggests that even for k = M5, one has A5 ~ 10k, which 
is a large enough separation for the 5D EFT to be useful. Hence, whenever we need to be 
concrete, we will take k = M5 when discussing the strong warping limit, although cases 
where k is somewhat smaller can be considered. 

We can further constrain the parameter space based on phenomenological considerations. 
First, one must reproduce the 4D Planck scale, Mp ^ 2.4 x 10^* GeV. Second, as will be 
seen in the following sections, the scale of KK resonances is set by A; = ke~^^^\ We will 
assume that A; ~ 0(1 TeV) to open the window to test this scenarios at the LHC These two 
constraints determine M5 and one of the three dimensionless parameters mentioned above. 
One can exchange one of the two remaining parameters for the radion decay constant A,, (to 



be discussed in Section 3.1 ), plus a free dimensionless parameter that we will choose as k/m. 
We elaborate next on the first two relations, before discussing the properties of the radion 
field. 

2.3.1 The 4D Planck Scale 

We start by computing the effective 4D Planck scale from the underlying Lagrangian pa- 



rameters. For instance, in the model defined by Eq. (14) the 4D (reduced) Planck mass 
is 



e-2^(^)rfy = ^ J dx [cosh(v^x)] , (24) 

^To be more precise, the KK scale is set by the warped down curvature scale at the IR boundaries, 
^cff e~"^'-^\ where fceff = ^'(i) = s^fc, and sa is a model-dependent number smaller than one. See Eq. (28) 
below. 
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where the factor of 2 arises from the Z2 symmetry of the background about ?/ = 0. There 
are three important scales in the problem: the curvature fc, the width of the domain wall 
1/m, and the size of the extra dimension L. As will be shown in the next subsection, we will 
be most interested in the case where m <^ k, i.e. when the domain wall is thick in units of 
the curvature scale. In this case, we find that Mp is well approximated by 



2M| 
m 



-Km 



4V2k 



tanh 




X mL 



(25) 



From this expression we see that unless the extra dimension is so small that mL ^ ^Jrnjk^ 
the tanh is essentially equal to one and we have Mp f» i(n j^fi)^!'^ M| / ^km ~ lA9M^/y/km. 

For completeness, we also give an approximate expression that holds in the narrow domain 
wall approximation, m ^ k: 



Ml 



2M| 



2M| 



kL 



dz 



k 



2v^fc/' 



m— 1 



fm r- 

cosh — v2z 
V k 

'^k/m _ 



-^/2k/m 



-2kL 



In the limit that A;/m — )■ we recover the RS result Mp 



V2m_ 
<-^^^)Ml/k [T]. 



(26) 



2.3.2 Constraints from TeV Scale New Physics 

We now explain how to further constrain the model parameters. Imposing that the observed 



Mp be reproduced, and for given /c/Ms, Eq. (24) determines as a function of k/m and 
mL. Fixing also k = ke~^^^\ one finds the "required" warp factor from A(L) = In k/k. 
Comparing to Eq. (23 ) allows one to fix mL, for given A;/m. One can then find 0o/0o from 
Eq (g. 

At this point the model is determined by e = k/M^ and k/m. In the left panel of Fig. [l| 
we show lines of constant mL in the e—{k/m) plane, requiring that Mp = 2.4 x 10^^ GeV 
and k = 1 TeV. This shows how to reach the regions of small [k <^ M5) and large {k < M5) 
warping. In the right panel, we show what this implies for the ratio M^/Mp (recall that M5 
is regarded as the fundamental scale in the theory, while Mp appears as an effective scale 
characterizing the interactions of the zero-mode graviton). 

Of course, naturality considerations lead to the expectation that k and M5 should not 
be very different, and that k ~ M5 ~ Mp as in the Randall-Sundrum proposal. In Fig. [2] we 
plot A{L) and mL as a function of k/m, assuming that M5 = k. We see that A{L) remains 
essentially constant at A{L) ~ 34—35. The curvature scale varies between (0.5 — 1.4) x 
10^^ GeV in the shown range. At the same time 0o/0o varies from 0o/0o ~ 1-1 for k/m few. 



For m <^k, Eq. (25) shows that Mp becomes essentiaUy independent of mL once mL > yjm/k, while 



for TO ^ /c Eq. (26) shows that Mp becomes essentially independent of mL for mL > m/k; in cither of these 
limits fitting Mp effectively determines A; as a function of k/m. Then mL is fixed as a function of k/m from 
Eq. (l23l). 



9 



0.1 
0.001 

^ 10-"^ 
10-' 
10-' 
10-" 



1 ■ ■ 1 1 1 1 




y 


1 1 ■ 1 


- &=:lTeV 


mL = 0.6 / 




0.5 / \ 








04 ■ 








0.3 








0.2 



50 



100 
k/m 



150 



200 



1 

0.1 
a. 0.01 
~J 0.001 

lo-"* 

10-^ 

lo-*" 





1 1 1 1 ■■ 1 1 


1^1 1 1 


~k = \ TeV 


mL = 0.6 




















0.3 






__JI2_ 









50 



100 
k/m 



150 



200 



Figure 1: k/M^ (left panel) and M^/Mp (right panel) as a function of k/m for different 
values ofmL, and for fixed k = ke'"^^^^ = 1 TeV. 
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Figure 2: A{L) from Eq. 23) (left panel), and mL from Eq. (22) (right panel) as a function 
of k/m. We impose Eq. (24) assuming M5 = k and k = ke~^^ = 1 TeV. 



down to 0o/0o ~ 0.5 for A;/m ~ 200. The right panel of Fig. [2j shows that mL > 1 (we also 
point out that 35 ^ /cL ^ 100 in the range shown in the figures). However, large values of 
mL require (j)o/{\/2(f)Q) ^ 7r/4, so that in general we expect mL ~ 0{1), and therefore we 
need k/m ^ 1 to generate the hierarchy. This means that the physical space is expected to 
be restricted to the "central region" of the domain-wall solution. As discussed in the previous 
subsection, in this region the background solution is dominated by the quadratic term in the 
scalar potential, and can be simply modeled by the linear superpotential of Eq. 

The small warping regime can be obtained for any value of k/m with mL < 1, as seen 
in the left panel of Fig. [l] For k/m ~ 1, Eq. (23) implies that A(L) ~ so that 

there is virtually no warping. In this case, the Planck-weak scale hierarchy is directly set 
by k/Mp ~ {k/M^) x {M^/Mp) ~ 10-^^ as can be seen from Fig. fTl N For A;/m > 1 one 



^'^Note that this gives a scenario quaUtatively similar to the UED models, where 1/L ^ k ^ k ^ TeV, 
but A/5 ^ lO^'^ GeV. Since the couplings among KK modes are of the same order as those of the 0-modes 
(recall A{L) ^ 1), the 4D EFT breaks down at O(IOO) x k, as in UEDs. However, as in warped scenarios, 
one would expect that UV brane observables are under control up to scales of order M5. 
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Figure 3: Left 'panel: <\){y) in units of k^^'^ for the "sine" superpotential (solid, red line), and 
the approximation of Eq. |^7| j (dashed, blue line). Right panel: A'{y) in units of k for the 
"sine" superpotential (solid, red line), and the approximation of Eq. ( 28) (dashed, blue line). 



arrives at the interesting situation where the Planck-weak scale hierarchy arises in part from 
a sizable warping 1 < A{L) ^ (9(35) as well as from a hierarchy associated with a small 
curvature k <^ M5. 

2.3.3 The Limit of a Fat UV Brane: Benchmark Scenarios 

The previous considerations suggest that we focus on the regime k/m 3> 1, corresponding to 
a "fat UV brane". In this case, (f){y) and A{y) can be approximated by linear and quadratic 
functions of y as follows: 

0(2/)/0o ~ s^y/L, (27) 
A{y) ^ SAkyV{2L) , (28) 

and sa = tanh (^\/2mLj in the model defined by the 
]q. (|14|). In other models, and sa will have a different dependence 



tanh 



mL \ 
V2 



where = \/2tan ^ 
"sine" superpotential, 

on the microscopic parameters, but will still be numbers of order one, so that Eqs. (27) and 



(28) provide a general parameterization of "fat brane scenarios", which we will call "the 
linear regime". The accuracy of the linear regime for k/m ^ 1 can be seen in Fig. [3| where 
we plot (f){y) and A'{y) as a function of ky, for the exact expressions (15) and (16), as well 



as for the approximations (27) and (28). In the following we will make extensive use of 



the approximate expressions (27) and (28) to gain intuition (and relatively accurate simple 
expressions) for the spectrum and couplings of KK modes. 

Here we will summarize the parameters for two "benchmark scenarios" that we will use 
in some of the numerical studies that follow. Both are based on the "sine" superpotential, 
i.e. on the profiles of Eqs. (15) and (16). As we will see in the following, /cgfr = A'{L){-- 
plays an important role. 

• "Strong warping scenario": we take M5 = k, k = 2 TeV and k/m = 200, 



SAk) 



apart 



from the phenomenological constraints discussed above. In this case, one finds k 
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5.2 X 10^^ GeV, A{L) ^ 33.2, mL ^ 0.50, kesL ^ 62, 0o/0o ~ 0.47. For reference, the 
radion decay constant, to be defined in the next section, is = 4.4 TeV. 

• "Small warping scenario": we take M5 = 2 x 10^ k, k = 2 TeV and k/m = 200, 
apart from the phenomenological constraints discussed above. In this case, one finds 
A; ^ 1.8 X 10^ GeV, A{L) ^ 4.5, mL ^ 0.18, k^gL ^ 9, 0o/0o ~ 0.18. The radion decay 
constant is A,. = 2 x 10^^ GeV. This example hes approximately at the right end of 
the blue, lower curve of Fig. [T} 

Here we have taken a relatively large value for k/m = 200. This is only to make sure that 



we are in the linear approximation of Eqs. (27) and (28), since some of the simple analytic 



expressions to be derived in the following sections rely on being in the linear regime. However, 
it is possible to take a smaller hierarchy between k and m, and we provide formulas for the 
general case as well. We also point out that the most notable difference between the strong 
and small warping scenarios enters through the radion decay constant, A^,, and is therefore 
relevant for radion or KK-radion phenomenology. The specific "small" warping scenario is 
inspired by a scenario for KK-radion dark matter, to be presented in |29j . 

We should also add that there are further phenomenological constraints on the KK scale 
from EW precision tests, and possibly from flavor and CP violating observables. However, 
such bounds depend on the particular model under consideration. Given that our emphasis 
in this paper is on the general framework, we do not engage in the detailed study of such 
constraints. Nevertheless, we have chosen our "strong warping scenario" in such a way that 
the KK gluon mass, Mkk ~ 2.39fceff ~ 3 TeV (see caption to Table [I] at the end of this 
paper), is roughly consistent with the bounds studied in other warped scenarios. For the 
"small warping scenario" we allow a smaller Mkk ~ 2.19 kcs ~ 1 TeV (see caption to Tableji]) 
since we expect the bounds to be somewhat looser in that case (somewhat similar to the UED 
case). We stress, however, that the data we present in this work are largely independent 
of the precise value of Mkk, within a given gravitational background, and should be useful 
when more detailed studies of the bounds on the KK scale in given models are performed. 



3 Properties of the Radion KK Tower 

We now turn to the scalar excitations associated with the background described in the 
previous section, paying special attention to the lightest KK-parity even scalar (the "radion") 
and the lightest KK-parity odd scalar (that will turn out to be the LKP). 

3.1 Normalization and Radion Decay Constant 

We start by discussing the normalization of the radion/scalar KK modes in general. Follow- 
ing Ref. [16] we parameterize the 4D physical scalar fluctuations (after gauge fixing) in the 
5D metric/bulk scalar system by 

ds^ = e-2^(^)-2^("'^)77^,dx'^dx" - [1 + 2F{x, y)f dy^ , (29) 



12 



(30) 



where A{y) and (f){y) are the background profiles described in the previous section, and 

^{x,y) = ^^e'^(y^dy[e-'^^y^Fix,y)] . (31) 
In terms of the KK decomposition 

oo 

F{x,y) = J2U^)^n{y) , (32) 

n=0 

one can show [1Q\ that F.n{y) obeys the "KK radion-scalar" equation of motion 

(33) 



with boundary condition (in the hmit that the $ VEV is frozen at the boundaries) 



(F^ - 2A'F, 



n)\±L 







(34) 



guaranteeing the self-adjointness of Eq. (33). In order to see the orthogonality relations for 



the KK wavefunctions that follow from this, it is useful to remove the terms proportional to 
in Eq. (g by defining F^iy) = e'^^y^ F^{y) / (p' {y) so that: 



F''- 



A" + 3A" - 2A' 



. 



(35) 



The standard argument then leads to the orthogonality condition, j^^^dy e^^FmFn oc 6mn- If 
we use the relation 0'^ = 3M|yl", implied by the Einstein equations for the background, we 
get the orthonormality relation for Fn{y): 



^ Fm{y)F4y) 
.r A"{y) 



-2A{L) 



kesm. 



2 ^mn ) 



(36) 



where the normalization of the wavefunctions, including the presence of the KK radion 
masses m„, was chosen for later convenience. 

Replacing Eqs. (29)-(31) and the KK decomposition of Eq. (32) in the action of Eq. ([T]), 
one finds that the Einstein-Hilbert term results in a contribution to the kinetic terms of 
the radion KK-modes, f„, given by 3M| Jd^x e~'^^FmFndfj,rmd'^rn, where the derivatives 
are contracted with the Minkowski metric. Similarly, the scalar part of the action gives 



a contribution to the KK radion kinetic terms of J d 



^ 1 

X I e 



-2A 



fnfmd^rmd^rn, with 



^^Sometimes it is useful to write this equation as i/)'^ ^tt dy [e ^"^F„] | + (e^"^m^ — 2A") F„ = 0, which 
suggests the relation between the existence of a massless radion mode, with Fq oc e^^, and A" = 0. 
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Figure 4: A^/Mp as a function of k/m for several values of mL. 
10^^ GeV is the reduced Planck mass. 



Recall that Mp = 2.4 x 



3jV/| OA 



e^"^ dy(e "^^Fj/). After integration by parts with respect to and using the KK radion 



equation (33) together with the relation (^'^ = 3M|A", this generates a term that precisely 
cancels the Einstein-Hilbert contribution, leaving behind only the terms proportional to m^. 
The action Eq. Q at quadratic order in f„(x) then reads: 

3 



S 



; ' - ---J A"(y) 

Al jd'xY,\dM^Wfn{^) + - 



(37) 



where we used the orthonormality relation (36), while the dots stand for higher order or 



non-derivative terms, and contain the potential for r„. We also defined the radion decay 
constant 



A^ 



'3M| 



-A{L) 



teff 



In the case that Eq. (|25| holds, i.e. for m < A;, we have A^ ^ (f )^^^Mpe 



(38) 



-ML) 



As remarked in Sectioii |2^ the hierarchy between k and M5 affects mostly the value of A^, 
for fixed k = ke~^^^\ As e = k/M^ decreases. A,, increases, while for k ~ M5, one has 
A^ ~ k. Thus, we see that the couplings of the radion are very sensitive to the size of the 
5D curvature. In Fig. |4| we show Ar /Mp as a function oi k/m for several values of mL, in 
the "sine" model of Eqs. (22) and (23). This plot should be compared to Fig. [l| but it is 
important to note that A^/Mp is actually independent of k. 
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3.2 Radion Mode 



In this subsection, we study the hghtest excitations of Eq. (33). Since in the hmit that 
A"{y) = there is a zero-mode solution with wavefunction Fq oc e^^, it is convenient to 
write in the general case 



(39) 



so that the radion equation of motion and the boundary conditions reduce to 

. 2{A'-^]q' = 2A" - e'Vo , (40) 

. (41) 



9 +9 



9 



±L 



Note that the differential equation is of first order in g' . We can obtain a good approximation 
to the ultralight mode with mass uiq by neglecting the g'"^ term, in which case Eq. (40) is 
solved by 



9'{y) 



(yy 



y A" 
2 / dz^e^^ 
0'' 



dz 



(42) 



where we chose the constant of integration so that g'{Q) = (since the radion is even about 



y = 0). The boundary condition (41) then implies 



mo 



p2A{z) 



dz 

4,'(zy 



.g4A(z)-2^(L) 



Cefre 



-ML) 



(43) 



which is a rather general expression for the "radion" mass in terms of (j){y) and A{ii). Here 
we have expressed the result in terms of k^a = A'{L) = SAk. 

For instance, f or the "strong" ("small") warping benchmark scenarios defined at the 
end of Subsection 



2.3.3 



which agrees very well with the numerical solution obtained without any approximations 



we obtain from Eq. (43) that mo = 0.22/ceff (^o = 0.67fcefr) 



eff 



effe 



-A(L)^ 



The neglected g'"^ term contributes negligibly to the mass because 



(here A;, 

it is comparable to the other terms in Eq. (40) only for y's near the UV brane, but the 
wavefunction profile, Eq. (39), is exponentially localized near the boundaries ai y = ±L. 



Hence distortions near the UV region can only give an exponentially small effect. 

As explained in Section [2| in order to solve the Planck/weak scale hierarchy problem 
the size of the fifth dimension is often stabilized in a region where (piy) and A{y) are well 



approximated by linear and quadratic functions in y, as given in Eqs. (27) and (28), re- 



spectively. Note that the radion EOM, Eq. (33), is approximately independent of the scalar 
profile, since (f)" in this limit. In this case, we can write a simpler expression for mg that 
is independent of the scalar profile: 



mo 



,2A(z) 



«efre 



-A{L) 



(44) 
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effC 



^A{L) 



(45) 



where in the second hne we used that when Aiy) is given by Eq. (28), the ratio of integrals 
in Eq. (44) is very close to 2. For the "strong" ("small") warping benchmark scenarios 



of Subsection 
mo ~ 0.25 /ceff 



2.3.3 



(mo 



Eq. i\U 
0.67 hs 



gives mo = 0.2Qkes {rriQ = 0.68 /ccfr), while Eq. (45) gives 



For future reference, we also give an approximate expression for the radion wavefunction. 



normalized according to Eq. (36), in the limit that A{y) is well approximated by Eq. (28). 
The fact that the radion mass is parametrically smaller than the KK scale, k^s, corresponds 
to the fact that g{y) in Eq. (39) is a small perturbation. Therefore, setting g{y) ^ 0, and 



using Eq. (45) for mo, we obtain 



where we used J^dpe^^ 



1 7'' 



Fo{y) ^ e2[^(^)-^(^)l , 
which holds whenever z ^ 1. 



(46) 



3.3 Lightest Odd-Scalar 



As already mentioned, the radion profile, Eqs. (39) or (46), is exponentially localized near 
the boundaries at y = ±L. Thus, after normalization, its value at ?/ = is exponentially 
small (while -being even- its derivative vanishes at the origin). One can then see that by 
adjusting mo by an exponentially small amount one can obtain a solution that vanishes 
at the origin, and therefore that there is an odd mode with a mass exponentially close to 
the radion mass given in Eq. ( [43| . We can make the argument more concrete by writing 
Fodd{y) = F{y) + e{y), where F{y) is the radion wavefunction given in Eq. (39). Writing 
also m 



n=l 



m 



+ Sm , where mo is the radion mass given above, Eq. (33) becomes 



2e' (^A' + + e e'^ {ml + 6m') - AA" + AA'^ 



e'^Sm'F = 



(47) 



We can approximately solve this equation in the limit where (f){y) and A{y) are given by 



Eqs. (27) and (28), respectively, which corresponds to the limit we are interested in. In this 



case, the terms proportional to 0" vanish. In addition, using the approximate expression for 
mo given in Eq. (45) we see that m^ e'^ / (AA") ~ e-'^['^iL)-My)]^ which is much smaller than 
one except for y ^ L. Since, as we will see, also Sm'^ <^ m^, we can focus on solving 



2A'e' - 4A"e + e'^6m'F 







(4J 



with A given by Eq. (28). Its solution takes the somewhat cumbersome form 
e{y) ^ e'^'^y^dm' ly / dz2F{z2 



Jy L J Z2 



-2A{2i) 



^^If the linear approximation is not excellent, one can obtain a pretty good approximation to the radion 

fc(,g e~^^^\ which only assumes that, due to the strong localization of the inverse 



mass by using mp f 
warp factor near y 



L over a distance of order 1/fcoff, the ratio of integrals in Eq. (44) is of order one. 
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(49) 



where we imposed e(L) = t'{L) = 0. The first condition can be obtained by requiring that 
-^odd(-^) = F{L) by a simple overall rescaling of the odd wavefunction, and the second follows 
from Eq. (34), together with the fact that the radion wavefunction, F{y), already satisfies 



this boundary condition. We emphasize that Eq. (49) assumes that A(y) is given by Eq. (28) 
Requiring that Fodd(O) = -F(O) + e(0) = fixes as 



6m 



d7 7p ^2\A{z)-A{L)] F{z) \ 



h^L F{L) "^^^ 



-2A{L) 



(50) 
(51) 



where we used ^(0) = ^'(0) = to write Eq. (50), while Eq. (51) holds due to the strong 
localization of the radion wavefunction and the inverse warp factor near y = L over a distance 
of order l/fceff- This expression shows that 5m? /m^ is (exponentially) small due to the IR 
localization of the radion wavefunction, Fiy). Hence, the radion and the lightest odd-scalar 
are exponentially degenerate in this scenario. 



4 Fermions 

We analyze now the fermionic sector. The Hermitian fermion action is 

^fermion = j d'x^ | '-'^e^JV^ D - '-{D M'^)^r%'jT^^ - y*<l>^^| , (52) 

where = {•y^, —^75) are the fiat space Dirac gamma matrices in 5D space, is the 
fiinfbein. Dm is the covariant derivative with respect to the gauge symmetry as well as 
general coordinate and local Lorentz transformations and yq, is a Yukawa coupling (with 



mass dimension —1/2). In the scalar background {(f)) given in Eq. (15) the Yukawa interaction 
gives rise to a y-dependent mass term for the fermion, mo = 

Note that the Yukawa interaction above is consistent with the Z2 transformation $(y) — 
— $(— y) provided either '^liy) ^hi—y) and '^niy) — )• —'^R^—y) or viceversa. This 
matches with the fact that, at a given KK level, the y — )■ —y transformation properties of 
the wavefunctions /£(?/) and /]^(|/) are opposite, but the associated 4D fields are nevertheless 
assigned the same KK parity (even for even n and odd for odd ra, independent of chirality). 
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For a diagonal metric of the form ds^ — a{yf'rj^^dx^^dx''' — h{y)'^dy^ the spin connection in Dm cancels 



out in the fermion action, Eq. ( 52 ) 
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4.1 Equations of Motion and Kaluza-Klein Decomposition 

Considering the action to quadratic order in the fields, the KK decomposition for each 4D 
chirahty of the 5D fermion reads 



n=0 



where the fermion profiles obey the orthonormality conditions 



(53) 



1 

2L 



dyflR{y)fTAy) 



(54) 



and the first order equations 



1 V 



The boundary conditions &i y = —L and y = L are given either by 



1 



fr+ nio--A'\l 



y=TL 



or 



fR-{mn + U']l 



, f^iTL) = 



, f2{TL) = 



y=TL 



(55) 
(56) 



(57) 



(58) 



where ' = dy. These give rise to a chiral zero-mode sector: Eq. (57) selects a left-handed (LH) 
zero-mode while Eq. (58) allows only a right-handed (RH) one. The zero-mode solutions are 
explicitly given by 



flniy) = Nexp(^^Aiy)T J'dzmniz)^ 



(59) 



where N is determined by Eq. (54). These solutions are even about y = since A{y) is even 
and the integral of the odd function m£,{z) is also even. 

We note that when m£,{y) and A'{y) are proportional to each other, which is the case 
when 0(?/) and A{y) are approximately given by Eqs. (27) and (28), the fermion zero-mode 
wavef unctions further simplify. In this case, we define a dimensionless parameter c by 



mniy) = ±cA'{y) 



(60) 



*We use the convention that Js^^l.r = T^L,fl- 
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where we take the upper (lower) sign when the 5D fermion has a LH (RH) zero mode (i.e. 
when the boundary conditions (57) or (58) hold, respectively). In our case, we have the 
explicit relation 



±l/*(s9i/SA)(0o/^) 



(61) 



where s^j, and sa were defined in Eqs. (27) and (28). In the limit that the proportionality of 



Eq. (60) holds, the zero-mode profiles of Eq. (59) take the simple form 



f{y) 



Noe 



(62) 



independently of the chirality of the zero-mode sector. This shows that for c ~ 1/2 the 
zero-mode profiles are (approximately) fiat, while for c > 1/2 (c < 1/2) they are localized 



near the UV brane (IR boundaries). Even in cases where Eq. (60) does not hold exactly, it is 



useful to formally exchange the Yukawa coupling for a c-parameter, defined by Eq. (61) 



and we shall do so from now on. The model-dependent constants, and sa can be easily 



estimated in the "central region" of the kink scalar profile [and were given after Eq. (28) 
for the "sine" superpotential, Eq. (14)]. In Fig. [sj we show the exact fermion zero-mode 
wavefunctions for several values of c to illustrate the situation in the benchmark scenarios 
defined in Subsection 2.3.3 , with k ^ m. This allows us to characterize the physics of 
localization in close analogy to the standard treatment of bulk fermions in the Randall- 
Sundrum scenario. But note that, in general, Eq. (59) does not admit an exactly fiat 
( "conformal" ) zero-mode solution, since the y-dependence of the two terms in the exponent 
cannot cancel for all ?/'s unless Eq. (60) holdsj^ Also, the UV-locahzed wavefunctions 
have a Gaussian-like profile (see also |l6l HZ]), while the IR-localized wavefunctions have an 
exponential-like profile a la RS. Nevertheless, it is interesting that we end up dynamically 
in a limit where the RS intuition qualitatively holds. The approximate expression for the 



zero-mode profiles, Eq. (62), will be very useful in the following sections. 



Going back to the heavy KK-mode fermion solutions, combining Eqs. (55) and (56) we 
obtain the decoupled system of second order differential equations: 



f2" - 2A'f2' + 
- 2A'f-' + 



^A'' - ^A" - moA' + m'j, - ml + e^V^ 



4 



-A 
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' All , A I I 2 I 2A 2 

-A + ttldA —nij^ — mj^ + e m„ 



fn 

Jl 

rn 
JR 



0, 



. 



(63) 
(64) 



In the absence of closed solutions for our chosen A{y) and moiy), one can find the spectrum 
(and wavefunctions) numerically by the "shooting method 
on 



For instance, concentrating 

j^, one can integrate numerically Eq. (63) by starting at y = —L with f2{—L) = 1 and 
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its derivative given by the first equation in (57). Varying x„ = nin/lkcse 



-ML)] 



one can 



^^This also shows that our background would break supersymmetry in a "hard way" 
zero-mode profile, hence its couplings, cannot exactly match those of the gauge field. 
^^From a numerical perspective, it is easier to solve for e^s^d') /£ ^j. 



since the gaugino 
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Figure 5: Fermion zero-mode wavef unctions, Eq. (59), for the "strong warping" (left panel) 
and "small warping" (right panel) scenarios. The c-parameter is defined by Eq. (61). Note 
that the Randall- Sundrum intuition approximately applies, but there is no symmetry between 
UV and IR localized wavefunctions due to the approximately Gaussian warp factor. The 
dynamical UV brane is at y = 0, while y = ±L correspond to IR boundaries. 



find those values of x„ for wfiich /]^'(0) = 0. These correspond to the even mode solutions. 
Similarly, the values of x„ that make /£(0) = determine the odd mode solutions. In 
practice, it is easier to find the even LH spectrum by solving Eq. (64) for f^ and requiring 
— which by the first order Eq. (55) implies /]^'(0) = 0. This only misses the LH 
zero-mode solution, that we know explicitly from Eq. (59). We show in Fig.lolthe result of the 



above procedure for the "strong warping" benchmark scenario defined in Subsection 2.3.3 



taking c = 1/2, so that the fermion zero mode is (approximately) flat. We use the exact 
profiles for A{y) and (j){y) to generate the figure. However, to the extent that Eqs. (27) and 
(28) hold, this example corresponds to the "conformal" case, and the solutions for the Xn 
approximately coincide with those of (unbroken) gauge fields]^ Indeed, one can check that 
the gauge spectrum coincides with the one shown in the figure at the percent level. We note 
from the figure the high degree of degeneracy between the even and odd KK modes. This is 
due to the fact that the KK-mode wavefunctions are highly peaked near the IR boundaries, 
hence small changes in x„ allows one to change from /]j'(0) = to (0) = 0, turning an even 
mode into an odd mode (the same effect discussed for the radion/odd-scalar in Section [3]). 
The degree of degeneracy increases as A{L) increases. 

As discussed above, there is always an even zero-mode. In addition, in certain regions of 
parameter space, there are light vector-like odd-modes that can be treated analytically, as 
we discuss in the next subsection. 



4.2 Light Vector-like Odd Modes 

For concreteness, in the following we concentrate on the case of a left-handed zero mode, 
i.e. we impose the boundary conditions of Eq. (57) [our results for the light mass and its 



Assuming vanishing brane-localized kinetic terms |48j . 
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Figure 6: We plot /i?(0) (solid, red curve) and /l(0) (dashed, green curve) as a function of 
Xn = rnn/[kes e^^^^^], for the boundary conditions of Eq. [51] [which select a LH zero-mode]. 
The zeros of fnlO) (which imply zeros o//^(0)j correspond to the even mode LH solutions, 
while those o//l(0) correspond to the odd mode LH solutions. We choose c = 1/2 [approxi- 
mately flat fermion zero-mode in the linear regime of Eqs. (21) and (28)], so that the roots 
also approximately correspond to those of gauge fields (see text). In this example, the roots are 
at x^™"^ ^ {0,2.43,5.51,8.60,11.70,14.80,...} and x"^"^ ^ {2.38,5.46,8.55,11.65,14.75,...}. 
Notice the near degeneracy between the even and odd spectra. 



wavefunction also hold when the zero-mode is right-handed, by a simple exchange of the 
labels L R]. In the case that the zero-mode fermion is exponentially localized near the 
IR boundaries, we expect to find an ultralight odd-fermion solution (for the same reasons as 
for the radion case discussed above). To find this small mass we write the wavefunction for 
the LH chirality of this odd-fermion (the first KK-mode), on the interval [0,L], as fl{y) = 
fiiy) + ^iy)y where /° is the (even) zero-mode profile given in Eq. (59) [fl{y) is completely 
fixed on the interval [—L, 0] from the fact that it is odd about y = 0]. One then finds that 



Eq. (63) becomes 



2A'e' + 



A'^ A" - rrioA' + m'^ - mj^ + e m' 



e + e m 



iJ L 



, (65) 



and requiring that the zero mode and the first excited mode have the same value at L we 
have the boundary conditions e(L) = e'(L) = 0. 

We give an approximate solution for e{y) in the limit that (f){y) and A{y) are approxi- 



mately given by Eqs. (27) and (28), respectively, and therefore Eq. (60) approximately holds 
(with the upper sign). Neglecting the e^'^m^e term (which is smaller than the other terms 
proportional to e provided mi <^ ke~^^^'^), Eq. (65) can be written as 



a 



-mlfl 



(66) 
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This can be integrated immediately to give 



Z2 



where we imposed the boundary conditions given after Eq. (65). Requiring that /^(O) 
/l(0) + ^(0) = determines the hghtest odd-fermion mass as 



where we used ^(0) = 0. This form shows that the hghtness of this odd mode is associated 
with /{J(0)//°(L) ^ 1. Indeed, since in the present hmit the zero-mode profile is well 



approximated by Eq. ( 62 ) , we can write 



mi 



This expression allows us to understand the suppression in mi compared to the natural KK 
scale fceff e"'^^'^-' when c < 1/2, i.e. when the zero-mode is IR localized and we expect a 
light odd-fermion mode [thus validating the approximation made above of neglecting the 



e m\e term in Eq. (66)]. In such a case, and using the approximation (28) for A(?/), the 
factor g- (2c- 1)^(21 ) jg localized near the IR boundary over a distance of order [(1 — 2c)A;efr]~^; 
and therefore the zi integration can be restricted to the region [L — {(1 — 2c)kcs}~^, L\. To 
proceed further, it is useful to distinguish two cases: 

• — 1/2 < c < 1/2: the factor e{2c+i)A(^2) jg localized towards the IR boundary, over a 
distance of order [(2c + l)A;efr]^^. We can therefore restrict the integral over Z2 to the region 
[L — {(H-2c)A;eff}~^5 L\. With this restriction over Z2 we find that the effective support of the 
Zi integration is of order (2fcefr)~^; in the vicinity of zi = L. Evaluating all the exponentials 
at L, we then get 

mi ^ V2(l + 2c)fceffe-^(^) , (70) 
where we recall that k^s = A'{L) = SAk. 

• c < —1/2: the factor e(2c+i)A(22) jg j-^q^ localized near ?/ = 0. The integral over zi can then 
be evaluated from to L for all Z2^ giving a factor [(1 — 2c)A;cfr]~^e~^'^'^*-^\ For |c| ~ 0(1), 
the integral over Z2 can then be estimated as an average and gives Le~"^*^'^''(l/2)[e^^'^"'"^'''^^'^'' + 
g(2c+i)A(o)j_ Putting these together, we have 




where we used ^(0) = 0. This shows that for extreme localization of the zero-mode solution, 
i.e. for 2c + 1 very negative, the lightest odd-fermion mode has an exponentially small mass, 
where mi/(A;eff 6"^^^^) ~ e(''+2)^(^). 
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Figure 7: Lightest odd-fermion mass ( in units of keg e'"^^^^ ) as a function of c. The solid 
(red) line corresponds to the numerical result with the exact scalar and warp factor profiles 
of Eqs. (15) and (16). The solid (blue) line corresponds to the numerical solution but with 
(j){y) and A{y) given by Eqs. (21) and (28) (the "linear regime"). The dashed (brown) 
line corresponds to the approximation of Eg. ("Approx. 1"). The dotted (green) line 



corresponds to Eqs. (70) and (71) ("Approx. 2"). 



In Fig. I?! we show the hghtest odd-fermion mass (in units of k^g e ^(^)) as a function of 
c for the approximation given in Eq. (69) with A{y) = kggy^ /{2L) [dashed, brown hne], and 



for Eqs. ( [70| and (71) in the regions —1/2 < c < 1/2 and c < —1/2, respectively [dotted, 
green hne]. ^® We also show the exact result (with the exact warp factor and scalar profile) 
obtained by numerical integration (solid, red line), as well as the exact numerical solution but 



assuming that (f){y) and A{y) are precisely given by Eqs. (27) and ([28|) (solid, blue line). We 



restrict to the "strong warping" scenario, with k/m = 200 (see Subsection 2.3.3). We see that 
the simple expressions derived above are fairly reasonable whenever mi/ {keg e~^^^^) % 1. 

For later use, we also give here the wavefunctions for the two chiralities of the first (KK- 
parity odd) fermion KK mode. The left-handed chirality (or more generally, the chirality 
that coincides with that of the zero-mode fermion in the tower) was found before in terms 



of e{y) of Eq. (67) on the interval [0,L]. Allowing for a normalization factor A^^i, to be 



determined from Eq. ( 54 ) , we have for y > 



fliy) ^ iVie(i-'=)^(-^) 



r(o) 



dzi 
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/°(0) 



(72) 



where f^ is the zero-mode wavefunction as given in Eq. (62). For y < one uses fl{—y) = 
—fliy). The wavefunction for the right-handed chirality (more generally, the chirality oppo- 
site to that of the zero-mode in the tower) is obtained from Eq. (55). Using that e^2~'=)^(y) 



^®The approximation of Eq. (70) fails at c = —1/2, since the factor e(^'^+^)"*(^2) becomes constant. However, 
numerically the approximation works rather well down to c ~ —0.48. For the dotted green line in Fig. [7] we 



interpolated linearly between Eq. (70 1 at c = —0.48 and Eq. (71) at c = —0.5 
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is the massless solution of Eq. (55) [in the hnear hmit], is given on the interval [0, L] by 

a+c)A(y) f^j.. f /°(^) 



/A(y)^-iVimie 



dz 



/°(0) 



(73) 



which exphcitly satisfies fji{L) = 0. For y < one uses ffj{—y) = +fR{y), i-e. that this 
wavefunction is even under y — )■ —y. To the extent that mi <^ kcse'~^^^\ we have that 
e{y) is a small effect and therefore the normalization constant A^i coincides approximately 
with the normalization of the zero mode: 



Nn 



2 (1 - 2c) A{L) 

g(l-2c)A(L) _ I 



(74) 



where we used the approximate expression for the warp factor, Eq. (28) [see Eq. (86) below 
for more details of the evaluation of A^ol- 



5 The Higgs Field and EWSB 

The action for a (for simplicity, real) bulk scalar H takes the form: 

^-scalar = j X ^ S^^Q^''' Om H On H - V (H) + 5iy + L) L + 5iy ~ L) C+L^ ,(75) 



where we allow for IR brane localized terms C±l (these (5-terms should be written with 
the induced metric, but since in our background we have ^/g = ^/gind, we have omitted this 
distinction for notational simplicity). Focusing on the quadratic terms, V{H) = ^M^if^+- ■ -, 
the KK decomposition is written as 



H{x^, y)= J2 Hnix^fniy) , (76) 
V 2-iv „ 

* n=0 

where the KK wavefunctions obey 

f';-2A'f',+ [A"-3A''-M' + e'^ml]U = 0, (77) 
and are normalized according to 



2L 



dyfn{y)fm{y) = 5nm ■ (75 



The boundary conditions associated with locahzed mass terms C±l = —\M±lH'^, where 
KK-parity imposes M_l = M+l, read 

/; + (A' ± M±^) /„ = 0. (79) 
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Since at y = we have m^i ~ and A' ~ 0, the fermion equations of motion (55) and (56 1 imply that if 
one of the two chiraHties has vanishing first derivative at y = (an even mode) then the opposite chirahty 
vanishes at y = (an odd mode), and viceversa. 
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It is convenient to parameterize the bulk and brane masses by 



{a^ - 4) k, 



2 

eff ) 



M. 



±L 



[a 



2) fccff + m±L ■ 



^0) 



where fceff = A\L), a = c + 1/2, and c = mci{L)/A'{L) as defined in Eq. (61) [see also 
Eqs. (27) and (28)]. In the AdSs limit, where A{y) = ky, these equations admit a zero-mode 



solution when m±L = 0. This follows by comparison with the KK fermion Eq. (63), and the 



boundary condition (57). The zero- mode localization is controlled by the same c-parameter 



that characterizes the fermion case. 

We point out that for modes that are highly localized near the IR boundaries, so that 
in the region y > we can replace A'{y) ^ A'{L) = k^s and A"{y)/k'^f^ ^ l/^k^sL) ^ 
1/[2^(L)] < 1, Eq. (fTTl) reduces to the fermion Eq. (g with mniy) ^ moiL) = cA'{L). 



The scalar and fermion boundary conditions also coincide in this case, up to the terms 
m±L. Thus, when c ^ 1/2 (strong fermion zero-mode localization in the IR), the scalar 
EOM also has a light mode, whose mass is controlled by m±L- To the extent that m±L is 
small, the scalar and fermion spectra almost coincide in this extreme IR localized limit. In 
particular, from the same considerations discussed in the radion and fermion cases of the 
previous sections, there is a KK-parity odd mode that is exponentially degenerate with the 
lightest scalar (KK-parity even) mode. 



5.1 Electroweak Symmetry Breaking 

A natural solution to the hierarchy problem suggests that the Higgs field is localized near the 
IR boundaries. In this case, the Higgs KK modes become heavy and decouple, except for the 
lightest KK-parity even and odd modes. If the Higgs had a vanishing VEV, the discussion of 
the previous section would imply that these two light states are nearly degenerate, with their 
masses controlled by the boundary mass terms m±L introduced above. However, we need to 
introduce a potential that leads to EWSB. It must also not break KK-parity, which implies 
that only the KK-parity even Higgs mode can acquire a non-zero VEV. The potential can 
arise from bulk or from IR localized terms, but we may simply analyze the physics of EWSB 
in the low-energy 4D effective theory. The relevant 4D degrees of freedom are KK-parity 
even and odd SU{2) doublets with hypercharge -|-1, that can be parameterized, respectively, 
as 

= ^.^^ + TJG" ) ■ { A- ) • 

where v denotes the warped-down VEV. Besides the kinetic terms, D^H\,D^H^+D^H^_D^H_ 
we can get potential terms of the form 

V = m^,,i7li/_ + ^Ai(i7li7_)2 + ^A2(i/ii/+-^)2)^ + A3(i/li/-)(i/li/+-^)^) 



(82) 



25 



where rn^^^ is the hght KK mass discussed in the previous subsection. The above potential 



terms lead to masses m^^ 



2A2'5 , mi 



m 



odd 



+ (A4 + As)^ 



"2, ml 



m 



odd 



+ (A4 - A5)f 



,~,2 



and mjj± 



m 



odd 



{G^ and are the eaten would-be Goldstone bosons). 



We note here, for later reference, that the KK-parity odd "kinetic terms" , D^H\ 
D^,H-D^'H+, as well as the KK-parity odd potential terms, \q{hIh^ + hIh+){hIh+ + 
(^H}_H^ — i)^) can also appear, for instance when coupled to a KK-parity odd radion excita- 
tion, r'/ A,,.. For strongly IR localized Higgs fields we expect Ai ~ ■ ■ ■ ~ Ag and C ~ 1- 



5.2 Yukawa Couplings to the Higgs and SM Fermion Masses 

We consider next the fermion Higgs interactions. Consider the top Yukawa interaction 



6S 



Y^nHQt + h.c. 



(83) 



where Y^d is the 5D Yukawa coupling (with mass dimension —1/2), H is the bulk Higgs 
doublet, Q is the third generation quark SU{2) doublet, and t is the top SU{2) singlet [Q 



satisfies the b.c.'s given in Eq. (57) and t those given in Eq. (58)]. 

Focusing on the zero-mode fermions, as defined by Eq. (53), the effective 4D Yukawa 
coupling reads 



Y, 



5D 



4D 



(2L)3/2 



dyhAy)fLiy)fRiy) 



5D rO 



fl{L)fn{L) 



where we took into account that the fermion zero-modes, as well as fhj^-, are even about 
?/ = 0, and the second equality holds in the limit where the lightest Higgs state is strongly 
IR localized over a distance of order fegfr, so that fh+{.L) ^ y/k^sL. We also defined a 5D 
Yukawa coupling Fs^) = Y^j:,/ y/2k^, as would be appropriate for a localized Higgs field. The 



general fermion zero-mode profiles are given in Eq. (59), but here we use the approximations 



given in Eq. (62) with the warp factor given in Eq. (28). Recall that we use conventions 



such that for c > 1/2 (c < 1/2) the zero- mode fermion is UV (IR) localized, independently 
of chirality. In addition, it is useful to write simple analytical expressions for the properly 
normalized wavefunctions that hold in the limit that A {L) 1. For c > 1/2 we use the fact 
that for 1 one has J^dpe'^ ~ AiA"/2, so that Eq. (62) reads 



fiy) 



8{c-l)A{L) 



TT 



1/4 



-(c-iM(2/) 



for c > 1/2 



^5) 



For c < 1/2 one has instead J^dpe^^ ~ e^^ whenever z 3> 1, and therefore 



fiy) 



^4 il - c) A{L) 

p(l-2c)A{L) 



,(|-c)A(y) 



for c < 1/2 



(86) 
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Figure 8: Left panel: Four- dimensional Yukawa coupling as a function of cr, for fixed cl = 
0.54, k/m = 200 and A{L) = 33.2. The dashed (green) line shows the approximate value of 
the top Yukawa coupling which occurs for cr ~ —0.45. Right panel: cl as a function of cr 
for the observed top mass ~ Ij, assuming = N(./{Ati), with = 3. 



These approximations work extremely well except very close to c = 1/2, where the wave- 
function is (approximately) flat: f^{y) ~ 1. Then, for instance, if the LH top has a fiat 
profile {cl = 1/2), while the RH top is localized near the IR boundaries {cr < 1/2) the 4D 
top Yukawa coupling reads 

F4D ~ a/(1-2c^) k,sL X ^^ , (87) 



where we used Eq. (28) for A{L). 

We estimate the maximum value of Yr,D from Naive Dimensional Analysis (NDA) in extra 
dimensions with singularities |15]. Writing the 5D Yukawa coupling as Y^u = w/A^, where 
A5 is the cutoff scale of the 5D theory, NDA gives i^nda ~ h/Vh, where I5 = 247r^ and 
U = 167r^ are the 5D and 4D loop factors, respectively. Also, taking A5 as the scale where 
5D QCD gets strong gives 

A5L ~ h/N, , (88) 
where A^c = 3 is the number of colors. Therefore, we estimate 

- - (89) 



In the left panel of Fig. [8] we show the top Yukawa coupling as a function of cr for 
Cl = 0.54, taking Y^d as given in Eq. ([89|) [and in the "strong warping" benchmark scenario 



of Subsection 2.3.3 . We use the exactly normalized zero-mode fermion profiles, although the 
approximate formulas given above lead to very accurate results except in the close vicinity of 
c = 1/2. Notice that we can generate a 4D Yukawa coupling of order one with a fiat -or even 
slightly localized on the UV brane {cl ~ 0.54)- left-handed zero mode, and a right-handed 
zero mode near the IR brane {cr ~ —0.45). For these values. Fig. [T] shows that there is a 
KK-parity odd resonance of the SU{2) singlet top with a mass mi ~ 0.44A;eff. 
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In the right panel of Fig. [8} we show the required function of c/j, that reproduces 

the observed 4D top Yukawa couphng, again assuming that the 5D Yukawa couphng takes 
the NDA value of Eq. (89). We see that the SU{2) top doublet has ~ 0.5, except when 
cr approaches 1/2. 



6 Gauge Fields 

For completeness, we summarize here the KK decomposition for gauge fields (in the absence 
of brane kinetic terms) . The results apply for general diagonal metrics of the form ^ . The 
gauge action reads 

S = jd^'x^l^- \g''''g''''FMNFKL } + S^,, . (90) 

The gauge fixing term in the background of Eq. ^ is chosen as [19] 

S,,. = d'x {d,V^ -idy [e-'^V,] Y , (91) 

where ^ is the gauge fixing parameter, and the 4D indices are contracted with rj^y. The KK 
mode expansions then read 

V,,,{x,y) = -^Y.KA^)fv,M ^ (92) 

* n=0 

where 

dy[e-'^dJ-]+mlf- = 0, (93) 

and = dyfy/rrin- For SM gauge fields, the boundary conditions are taken as 5j^/y|o,L = 
f^\o,L = 0, so that there is no zero-mode for V5. In the above i?^-type of gauge, has mass 
^m^. The KK wavefunctions obey the orthonormality relation 



1 

2L 



"^dyf^f^ = 5^m. (94) 

L 



The gauge eigenvalues in the "strong warping" benchmark scenario defined in Subsec- 
tion 2.3.3| are approximately given in Fig. [6j The gauge zero-mode wavefunction, however, is 
always exactly flat {fy = 1), unlike the zero-mode wavefunction for a fermion with c = 1/2. 



7 Radion and KK-Radion Couplings at Linear Order 

In previous sections we saw that, apart from the "zero-mode" sector (that includes the SM 
fields plus the radion), some first-level KK states can be relatively light compared to the KK 
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scale. These light states are KK-parity odd. In particular, in the present scenario the lightest 
KK-parity odd particle (the LKP) is the first level KK-radion while the second lightest KK 
state (the NLKP) corresponds to either the first KK-parity odd mode of the IR localized top 
SU (2) singlet, or the first KK-parity odd excitation of the Higgs doublet. In this section, we 
work out the Feynman rules for the interactions involving a radion/KK-radion with the top 
KK tower, and give some numerical examples. We also give the radion/KK-radion couplings 
to the gauge KK towers, and also to the two Higgs doublets. 



7.1 KK-Radion Couplings to Fermions 



The couplings of the radion KK tower to the KK fermions are obtained by replacing into 



the action the KK decompositions for the metric and bulk scalar, as given in Eqs. (29)-(32) 



together with the fermion KK decomposition given in Eq. (53). We are interested in the 
terms linear in the radion or its first KK-parity odd excitation. "^^ There are two types of 
contributions. First, those arising from the bulk fermion action, Eq. (52), which involve KK 
fermions belonging to the same 5D fermion field. Second, those arising from the localized 
Yukawa interactions, as in Eq. (83), which involve KK fermions belonging to different 5D 
fermion fields. 

The Yukawa interactions lead to mass mixing (after EWSB) between different fermion KK 
towers, most importantly for the top tower. Although these effects can be taken into account 
perturbatively, it is safer to perform a KK decomposition that already incorporates EWSB 
effects (for an IR-localized Higgs field, the bulk solutions worked out before do not change, 
but the boundary conditions are deformed by EWSB). In this way, the effects of EWSB 
appear directly in the vertices and spectrum, and one does not have to worry about mixing 
in external legs of physical processes. A second possibility is to diagonalize the fermion 
mass matrix including a sufficiently large number of KK modes. Here, for simplicity, we will 
derive the Feynman rules for the vertices ignoring EWSB. This should contain the essential 
physics for most applications (but see Ref. [29] for a discussion of EWSB effects). However, 
note that the localized Yukawa terms also lead to quartic interactions involving the physical 
Higgs, h, of the form h—Ti—P—t^ that are proportional to the EWSB zero-mode fermion 
mass, rrit. 

As just stated, we proceed now by setting EWSB to zero. In terms of the canonically 
normalized radion KK modes, r„(x) = Arf„(x), where A^ is the radion decay constant defined 



in Eq. (38), the interactions induced by the bulk action take the form 



/'bulk 



oo 



ijk=0 



+ 



ijk=0 



\mf,^'n^i + h.c.] ,(95) 



^"^Note that the coupUngs involving heavier radion KK states receive a further suppression from the l/m^ 
in the radion normaUzation, Eq. (36) [in other words, heavier radion modes have an effectively larger decay 



constant, A" ^ Xn-^r, where m„ ~ Xn fccff e is the n-th KK radion mass]. 
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where iIj^x=\ [^l^d^x - ('9/,^)7''x] ■ 



LL 
{Jijk 



2L 



m 



RL 

ijk 



L 



RR 

yijk 



1 

2L 



(96) 



2L 



dy 



m, 



2mDFi + ^e'^(e- 



2 A 



FY 



fyi , (97) 



and m^^f 



m. 



RL 
ikj ■ 



Here we have used the fact that the funfbein is diagonal with 



e {1 + F) 5^ and ~ (1 — 2F). To arrive at Eq. (97) we used the fermion equations of 
motion, Eq. (55) and (56). The second term in the integral in Eq. (97) arises from the 
scalar fluctuations given in Eq. (31), where we used the background relation (^'^ = 3M^A" 
to eliminate M5. Note also that the integrands in Eq. (96) and Eq. (97) are odd when 

^LL ^RR 

'ijk Jijk 



{i + j + k) 7^ mod 2, which implies g^^ — 



m 



RL 
ijk 



in such cases, so that KK-parity 



is not violated. 

It is of some interest to consider the interactions involving n radion modes and a fermion 
pair. The ■ ■ -Ti^ip^ip^ interactions take a similar form to Eq. (95) with — )■ A", gf^^ — )■ 
Cngft.ir.jk, gfjk Cngf,^..i,,jk and mf^ dnmf^ .^.j^, where c„ ^-3)" (2n/3 - 1) /n! and 
dn = (— 4)"'~^/?T,! are combinatorial factors, and 



LL 

Htl...lnjk 



1 

2L 



dy F,, . . . FUift 



-L 



RR 

Htl...lnjk 



2L 



dyF^.-.-FUykA^^) 



m. 



RL 

ii...injk 



'^g^.^jk + 2gf^\^, mk-j{2~n) I dy e-^mnF,, ■ ■ ■ F^kfl 



fj fk 



1 n 
2L2 



(3 - n) £^ dy ^e^F, • • • {e-'^Fj'^f, . 



(99) 



It is understood that in the Lagrangian we write separate terms for each (zi, . . . , in) ordered 
ra-tuple. For instance, in the case n = 2, if zi 7^ 12 we have two terms corresponding to 
Ti^Ti^ip^ip^ and Ti^Ti^ip^il)^ , with the coefficients given above. If ii = 12 then there is a single 
term. This case plays a role in the annihilation processes TiTi — )■ ipil). 

The bulk action also leads to interactions involving the radion, fermion pairs and gauge 
bosons of the form 



Ijlk 



(100) 



where we used the KK expansion for the gauge fields, Eq. (92), and defined 

-L 



LL 

gijik 



RR 

gijik 



gv_ 

2L 

gv_ 

2L 



dyFJifyfl 



-L 
L 

( 

~L 



dyF^y^R 



(101) 
(102) 
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with gv = g5v/V'2L the 4D gauge couphng. We work in unitary gauge so that there is no 
i^L^b'^R interaction (the V5 components of the gauge field are eaten by the higher KK modes, 
and the zero-mode of V5 is projected out by the boundary conditions). Note that when we 
consider the interactions involving a zero-mode gauge boson, with fyi^y) = 1 up to EWSB 



effects, we have ^^^^4 = 9v9ijk gfj^ = gvQ^k^ where gf^^. and g^'^ are defined in Eq. (|96|) 

The Feynman rules for a single radion state, r^, and two KK fermions or a pair of 
fermions and a (neutral) gauge boson are then 



k 




s 




b 



A, 



which, as mentioned above, are non-vanishing only when {i + j + k) = Q mod 2, or {i + 
j + k + I) =0 mod 2. For the Feynman rule involving the gauge boson: in the case of the 
photon, {T"-)bc 1 and gy = Qt^, with = | the top quark electric charge; for the Z 
gauge boson, (T")^^ -> 1 and gy = \f(F+lf' {Tf - QtSw)] for the gluons, gy = y/^na~s, 
with T" the SU{3) generators normalized according to Tr(T"T^) = ^6""^, and a, b, c are color 
indices. 

The vertices involving the lightest states (the radion, the first KK radion, the RH top 
and its first KK excitation) are especially important, and their Feynman rules can be further 
simplified. We can also restrict to couplings involving the zero-mode gauge bosons. Since 
the radion, r [i.e. the lightest of the scalar fiuctuations in the 5D metric/bulk scalar system. 



with a mass mo as given in Eqs. (43)-(45)], is KK-parity even, it can only couple to a pair 
of KK-parity even or a pair of KK-parity odd fermions. Similarly, the first radion excitation 
[which we call r' and is essentially degenerate with r (see Section [3])], is KK-parity odd and 
can couple to t and t', where we denote by t' the first KK excitation of the SU(2) singlet 
top. The couplings involving these light fields can be approximated based on the following 
observations: 

1. The radion wavefunction is highly peaked near the IR boundaries over a distance of 



order l/k^s, as given in Eq. (46). Also, as discussed in Section 4.2, r' has a wavefunction 



essentially identical to that of the radion: Fi{y) ^ FQ{y). 

2. For the RH top, /^(y) ~ fniy) on the interval [0,L], and both the zero-mode and its 
first KK excitation, t', are highly peaked near the IR boundaries. Recall that f^{y) 



was given in Eq. (86). 
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3. The wavefunction for the LH chirahty of t', fl{y) [given in Eq. (73), with the trivial 
change in label i? — )■ L] is more spread-out over the full extent of the extra-dimension. 
To see this, note that the IR localization of f^{y) implies that the integral in Eq. (73) 
is almost y-in dependent, except very close to y = L, where it vanishes. Hence the main 
source of y-dependence arises from the e'-a^'^^^'^*-^-' factor. For cr ~ —1/2 this factor is 
also nearly constant, so that the normalization implies that fl{y) ~ 1. 



It follows that for cji ~ —1/2 the couplings involving /^(y) or /^(y) are larger than 
those involving fl{y) by a factor of order f^{L) ~ 1^2(1 — 2cji)A{L) ~ \^kcsL, for each 
occurrence of the RH versus LH wavefunction, so that for instance g^Q^ ~ g^fl^/ {k^sL) in 



Eq. (96). When cr is not so close to —1/2 the suppression is less severe, but we still have 



9i(k ^ fl'ioi provided f^iy) is IR localized. The integral in Eq. (97) is similarly suppressed. 
In fact, the second term inside the integral goes like e'^g'Fof'^fl, where g', as given in 



Eq. (42), gives an additional suppression that makes this contribution arising from the scalar 



fluctuations completely negligible. When focusing on the couplings to the lightest top KK 
state, we can therefore neglect all the contributions that depend on fl{y)- For typical choices 
of parameters the errors thus induced are of order (10—20)%. 

The upshot is that the bulk interactions between the radion/LKP, the NLKP (with both 
RH and LH chiralities) and the top quark, as well as those involving a SM gauge boson, are 
controlled by 



fi'ooo 



9on 



9wi 



RR 
i/llO 



1 -2c 
3 - 2c 



.1/2 



(103) 



For the radion couplings to a fermion pair, the explicit Feynman rules are: 




(104) 



A. 



(105) 



The mass of the first KK resonance, mf, is given in Eq. (70). Similarly, the couplings 
involving r' are: 




9101 

Ar 



k+^f)PR-2mt>PL 



(106) 
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(107) 



where stands for any of the neutral SM gauge bosons. We did not write exphcitly the 
hermitian conjugate processes. The rules for the Vn—r—t—i and V^— r— t'— f vertices are 



analogous to Eq. (107), with g^^^ — )■ and g^^^ — )■ g^^, respectively. 

For illustration, we give an example of the relevant couplings for the radion [g^^^, g^^ 
and mQ^^/fceff e~^*-^^] and for the KK-radion [(^f^^? 9ijk ^^"^ "^jj^/^efr e""^*-^^], to the lowest 
lying KK fermions. In Table [T] we give the couplings in the "strong warping" benchmark sce- 
nario, while Table [2] contains the corresponding couplings in the "small warping" benchmark 
scenario (these scenarios are defined in Subsection 2.3.3). We chose c = —0.2, as might be 



appropriate for the SU{2) top singlet. In this case, the fermion masses are given by 



/ p-A{L) 



with = 0,1.15,3.38,4.32,6.50,7.43, 



[X 



f 



0,0.97,2.96,3.78,5.63,6.51, 



in the strong (small) warping benchmark scenarios. Even (odd) n correspond to KK-parity 
even (odd) states. For such c-values, the KK-parity even and odd states are not as degen- 
erate as for c = 1/2, but one can recognize the corresponding semi-degenerate pairs. Note 
that the tables display the selection rules implied by KK-parity conservation. We also point 



out that the last integral in the definition of m^^ [see Eq. (97)] typically gives a small con- 



tribution, although it has been included in the numerical evaluation of Tables [T] and |2} We 
see in the tables that the couplings involving fermions other than the zero-mode and/or the 
first ("ultra-light") fermion excitation, are almost vector-like (see the 4x4 submatrices in 
the lower right corner of each table). The couplings involving t and t', on the other hand. 



are close to chiral, as we have also exhibited in Eqs. (106) and (107). 



7.2 KK-Radion Couplings to Gauge Bosons and the Higgs Fields 

We also summarize the Feynman rules for one or two radion KK modes and two gauge KK 
modes. We focus on the interactions with V^" which are sufficient for tree-level calculations 



(where one can work in unitary gauge with = 0). The F^yF^'-' terms in Eq. (90) give 

1 1 



4 2L 



Y^[l + 2F{x,y)]fyf'yFf,y{x)F^'^\x) 



(108) 



j,fc=0 

(with indices raised by the Minkowski metric), while the F^^F^^ terms give 



2 2L 



-2[A{y)+F{x,y)] 



[1 + 2F(x, y)]'' dyfydyf^ VUx)V^\x) 



(109) 



where we used the parameterization of the radion modes given in Eq. (29), as well as the 



gauge KK decomposition, Eq. (92). Using also the radion KK decomposition of Eq. (32) 



33 



with ri{x) = Arri{x), the Feynman rule for the vertex involving a single radion and two 
gauge fields reads 




-r,: 



2i 

A, 



(110) 



where 



dijk 



1 

2L 



j fk 
V 1 



n 

Uijk 



2L 



dye-'^FJif^, 



and we restored the gauge indices, a,b. Here we used fl = dyfy/rrij, since it is simpler to 
find fliy) numerically. Note that Eqs. (Ill) with Fj — )■ 1 are precisely the orthonormality 
relations for the fy and . For completeness, we also quote the interactions of two radions 



and two gauge bosons, which arise from Eq. (109), although these do not involve the gauge 
zero-modes: 




foft m 



(112) 



where g^^^ is defined as g^j^ in Eq. (Ill), but with two radion wavefunctions, FiFi. 

Finally, we give the Feynman rules for the Higgs interactions with one and two radion 
modes. Here we assume, for simplicity, that the Higgs fields are IR localized, so that the 
radion couplings enter through the induced metric. We also neglect the possible mixing 
between Higgses and radion modes. Referring to the discussion of kinetic and potential Higgs 



terms of Subsection 5.1, one must only remember that the radion modes enter through a 



multiplicative factor e "^^^^'y^ when coupled to scalar kinetic terms, and through e ^-^(^'2^) 



when coupled to potential terms. Using the radion KK decomposition, Eq. (32), and in 



terms of the Higgs mass eigenstates defined in Eq. (81 ), the Feynman rules involving a single 
radion mode are 



ni + (-ir] 



Ar 

A, 
F^iL) 



[ki-k2 + 2ml] , 



(113) 



(114) 



(115) 
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More concretely, in order to find the gauge KK spectrum, to,- , we apply the "shooting method" described 



at the end of Subsection |4.l| to = e 



'2A 



/| , which -from Eq. (|93|)- obeys +m'je 



2„2Ajy^0and fi 



= 0. 



±L 
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where = h^,h^,a,H'^,G^,G'^ (and 'm^QO q± = 0). We also denoted the (incoming) Higgs 
momenta by ki and k2-, and Fi{y) is the wavefunction. The Feynman rules involving two 
radion modes are 



+ [k,-k, + 4ml] , (116) 



r^r.h^h. -2. [1 - (-1)-+^] ^^^Y/^^ [kvk2 + Sm^] , (117) 

r 

r,r,{a,G')/iH^G^) — > -2t [l - {-ly^^] ^^^^^ k.-k^ . (118) 

8 Summary and Conclusions 

In this work we have introduced a 5-dimensional scenario where the dynamics that stabilizes 
the size of the fifth dimension: i) induces the non-trivial warping that allows understanding 
the weakness of the gravitational interactions compared to the weak interactions ii) leads to 
a metric background that implies the existence of a KK-parity symmetry, thus predicting a 
weak scale stable particle (the LKP), iii) leads to fermion localization via Yukawa interactions 
involving the stabilizing field, and iv) the Higgs sector is described by a THDM with an "in- 
ert" Higgs doublet (tan/3 = oo). It is argued that generically the LKP is the first excitation 
associated with the radion field. This KK-radion is expected to be highly degenerate with 
the radion, and its interactions are controlled by the same decay constant that controls the 
radion interactions. This allows to infer properties of the KK-radion from collider signatures 
of the radion, and help identify the collider KK-radion signatures. In particular, it may be 
possible to probe whether the LKP accounts for the DM content of the universe via collider 
measurements. Although we have not explored this question in this paper (we present the 
analysis in the companion paper [29]), we have presented the technical ingredients necessary 
to study the phenomenology of such scenarios. In particular, we provided simple analytical 
expressions for the masses and wavefunctions of the lowest lying KK states, which can be 
summarized as follows: 

• The KK scale is set by the warped down IR curvature scale, kes = k^s e'^^^'' = 
v4'(L)e-^(^), where 2L is the proper size of the extra dimension, while A{y) depends 
on the stabilization mechanism (and needs not be linear in y, as is the case of AdSs). 

• The radion (KK-radion) is strongly and symmetrically (antisymmetrically) localized 
near the IR boundaries, according to -Fo(l/) ~ e^'^*-^-'"^''^^' for y G [0,L]. Their masses 
are exponentially degenerate at tree level, and given approximately by k^s in 



the case that the dynamical "UV brane" is fat [see Eq. (43) for a completely general 
expression] . 

The zero-mode fermions have wavefunctions that behave like e"'''^"^)^^?') in the "fat 



UV brane" limit [see Eqs. (62) and (74), and Eq. (59) for the exact expression in a 
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general 4D Lorentz invariant background]. The effective c-parameter allows a simple 
characterization of the localization properties. In the "fat UV brane" limit, UV lo- 
calized wavefunctions (c > 1/2) are Gaussian-like, while IR localized wavefunctions 
(c < 1/2) are approximately exponential. 

The lightest KK fermion excitations that are localized towards the IR boundaries have 
a mass approximately given by mi ^ ^2{l + 2c) when -1/2 < c < 1/2. Their LH 



and RH wavefunctions were given in Eqs. (72) and (73) 



We also provided sample KK spectra and couplings between radion/KK-radion and 
fermion/gauge modes in Tables [l] and [2| and briefly discussed the fermion couplings to 
the Higgs, as well as the KK-radion/Higgs interactions. These ingredients are relevant when 
studying the phenomenology of the present scenario. We discussed the possibility that the 
curvature scale ranges from order the fundamental Planck scale down to much smaller values. 
The former case is similar in spirit to the RS proposal, while the latter has some features in 
common with Universal Extra Dimensions, but allows a consistent description of the gravi- 
tational interactions and, via the AdS/CFT correspondence may be a valid description up to 
a scale exponentially larger than the weak scale (but well below the 4D Planck mass). This 
latter scenario may allow for a lighter spectrum of KK states than is expected in RS models. 
The EW and flavor constraints deserve a more detailed study, as does the associated collider 
phenomenology of the Z2-symmetric warped scenarios. 
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Table 1: Radion (left) and KK-radion (right) couplings to KK fermion pairs [see Eq. (95) 
in the "strong warping scenario", for c = —0.2. We also show the KK-radion-gauge boson 
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boson ones are Xq 



couphngs [see Eq. (111)]. The lowest lying KK fermion masses, in units of kes = kes^ 
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0, x\ 



3.38, x{ ^ 4.32, xi ^ 6.50 and x 
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2.39, X 
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7.43, while the gauge 
2.41, xX 5.46 and xl ^ 5.49. The states 
with mass Xn with even (odd) n are KK-parity even (odd). In general, EWSB will lift Xq 
to a non-vanishing value. The radion and KK-radion masses are mo ~ mi ^ 0.22 k^s- The 
dimensionless values in this table were computed assuming k = 2 TeV (which corresponds 
to fceff ~ 1-2 TeV), but are rather insensitive to this choice, as long as ~ 0{TeV). 
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Table 2: Radion (left) and KK-radion (right) couplings to KK fermion pairs [see Eq. (95) 
in the "small warping scenario" , for c = —0.2. We also show the KK-radion-gauge boson 
couplings [see Eq. (111)]. The lowest lying KK fermion masses, in units of fcgff = kes^'^^^^K 



are Xq = 0, a;{ ~ 0.97, X2 ~ 2.96, ~ 3.78, 0:4 ~ 5.63 and ~ 6.51, while the gauge 
boson ones are x^ = 0, xj^ ~ 2.19, x^ ~ 2.41, x^ ~ 4.86 and x^ ~ 5.17. The states 
with mass x^ with even (odd) n are KK-parity even (odd). In general, EWSB will lift Xq 
to a non-vanishing value. The radion and KK-radion masses are mo ~ mi ^ 0.66 k^s- The 
dimensionless values in this table were computed assuming k = 2 TeV (which corresponds 
to keff ~ 500 GeV), but are not very sensitive to this choice, as long as A; ~ C(TeV). 
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